Cotunneling Transport and Quantum Phase Transitions in Coupled 
Josephson-Junction Chains with Charge Frustration 
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We investigate the quantum phase transitions in two capacitively coupled chains of ultra-small 
Josephson-junctions, where the particle-hole symmetry is broken by the gate voltage applied to each 
superconducting island. Near the maximal- frustration line, cotunneling of the particles along the 
two chains is shown to play a major role in the transport and to drive a quantum phase transition 
out of the charge-density wave insulator, as the Josephson-coupling energy is increased. We also 
argue briefly that slightly off the symmetry line, the universality class of the transition remains the 
same as that right on the line, being driven by the particle-hole pairs. 



PACS numbers; 74.50. -f-r, 67.40. Db, 73.23.Hk 

Systems of ultra-small tunnel junctions composed of 
metallic or superconducting electrodes have been the 
source of a great number of experimental and theoret- 
ical works Single charge (electron or Cooper pair) 
tunneling observed in those systems demonstrates the 
remarkable effects of Coulomb blockade. Especially, 
in Josephson-junction arrays, the charging energy in 
competition with the Josephson-coupling energy further 
brings about the noble effects of quantum fluctuations, 
which induce quantum phase transitions at zero temper- 
ature Very recently, another fascinating manifesta- 
tion of Coulomb blockade has been revealed in capac- 
itively coupled one-dimensional (ID) arrays of metallic 
tunnel junctions j^,^: In such coupled chains, the ma- 
jor transport along both chains occurs via cotunneling of 
the electron-hole pairs, which is a quantum mechanical 
process through an intermediate virtual state. Such a co- 
tunneling transport leads to the interesting phenomenon 
of the current mirror. 

In capacitively coupled Josephson-junction chains, the 
counterpart of the electron-hole pair is the particle-hole 
pair, i.e., the pair of an excess and a deficit in Cooper 
pairs across the two chains. Such particle-holes pairs, 
combined with the quantum fluctuations, have been pro- 
posed to drive the insulator-to-superconductor transition 
lH . Here it should be noticed that the particle-hole pair 
is stable only near the particle-hole symmetry line; far 
away from the symmetry line, it does not make the low- 
est charging-energy configuration anymore. Moreover, 
in a single chain of Josephson junctions, breaking the 
particle-hole symmetry (by applying a gate voltage) is 
known to change immediately the universality class of the 
transition . Therefore, it is necessary to find another 
relevant cotunneling process, if any, off the particle-hole 
symmetry line and to examine how the transitions change 
in coupled Josephson-junction chains. 

As an attempt toward that goal, we investigate in 
this paper the quantum phase transitions in two chains 



of ultra-small Josephson-junctions, coupled capacitively 
with each other. The particle-hole symmetry is broken by 
the gate voltage applied to each superconducting island; 
the resulting induced charge introduces frustration to the 
system. Near the maximal-frustration line, cotunneling 
of the particles along the two chains is found to play a 
major role in the transport and to drive a quantum phase 
transition out of the charge-density wave (CDW) insula- 
tor, as the Josephson-coupling energy is increased. We 
also argue that slightly off the symmetry line, the uni- 
versality class of the transition remains the same as that 
right on the line, i.e., a Berezinskii-Kosterlitz-Thouless 
(BKT) transition driven by the particle-hole pairs. 

We consider two ID arrays, i.e., chains of Joseph- 
son junctions, each of which is characterized by the 
Josephson coupling energy Ej and the charging ener- 



gies Eo = e /2Co and Ei 



/2Ci, associated with 



the self-capacitance Co and the junction capacitance Ci, 
respectively (see Fig. |l|). The two chains are coupled 
with each other via the capacitance C/, with which the 
electrostatic energy Ej = e^/2C/ is associated, while 
no Cooper-pair tunneling is allowed between the two 
chains j^]. The intra-chain capacitances are assumed to 
be so small {Ej <C Eo,Ei) that, without the coupling, 
each chain would be in the insulating phase jl^. We are 
interested in the limit where the coupling capacitance is 
sufficiently large compared with the intra-chain capaci- 
tances, Ci > Co,Ci, i.e., Ei < Eq,Ei [see Eq. (|) be- 
low] . On each superconducting island, external gate volt- 
age Vg is applied, and accordingly, the external charge 
Ug = C'oVg/2e is induced, with e being the electric charge. 
The external charge Ug breaks the particle-hole symme- 
try of the system, introducing charge frustration. We 
restrict our discussion to two regions: near the particle- 
hole symmetry line (|ng— A/"] <C 1/4 with JV integer) and 
near the maximal- frustration line (|ng— A/"— 1/2| ^ 1/4), 
where the properties of the system are severely different. 
We further note the invariance with respect to the sub- 
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stitution Ug — > rig + 1, and take = without loss of The Hamiltonian describing the system is given by 
generaHty. 



where the number ni{x) of the Cooper pairs and the 
phase 0f(x) of the superconducting order parameter 
at site X on the ^th chain (£ — 1,2) are quantum- 
mechanically conjugate variables: [ni{x), (j)e{x')] = 
iSu'Sxx'- The capacitance matrix C in Eq. (j^) can be 
written in the block form: 

Cu'{x,x') = C{x,x') J 5 +Sx,x'Ci _\ ~\ (2) 
with the intra-chain capacitance matrix 

C{x,x') = CqSxx' + Ci [2Sxx' - Sx,x' + 1 - Sx,x'-l] ■ 

For simplicity, we keep only the on-site and the nearest 
neighbor interactions between the charges (i.e., Ci/Cq < 
1) although this is not essential in the subsequent dis- 
cussion (as long as the interaction range is finite). With 
the block form of the capacitance matrix in Eq. (|^), the 
Hamiltonian can be conveniently expressed as the sum 

H = H"c+H^+Hj (3) 

with the components 

H'c = U„Y.[n+{x) ~ 2n,]2 + V„ Y.i^-{x)]' 

X X 

= UiJ2[n+ix) - 2ng][n+{x+l) - 2ng] 

X 

+ ViYn.{x)n.{x+l) (4) 

X 

l,x 

where n±(x) = ni{x) ± n2{x) and the coupling strengths 
are given by Uq ~ 2Ea, Ui ~ A(Ci/Ca)En, ~ Ei, and 
Vi ~ {Ci/Ci)Ei. 

The on-site charging energy term of the Hamiltonian 
i?^ in Eq. (^ reveals clearly the crucial difference be- 
tween the charge configurations in the system near the 
maximal- frustration line Ug = 1/2 and those near the 
particle-hole symmetry line = 0. In the former region 
(|ng— 1/2| <C 1), the charge configurations which do not 
satisfy the condition n+(x) = 1 (for all x) have a huge 
excitation gap of the order of Eq. (Note that we are in- 
terested in the parameter regime E'/, Ej <C Eq, Ei.) Fur- 
thermore, the ground states of H^, separated from the 



excited states by the gap of the order of Ej, have two- 
fold degeneracy for each x, corresponding to ri_ (x) = ±1. 
This degeneracy is lifted as the Josephson-coupling en- 
ergy Ej is turned on. As a result, it is convenient in 
this case to work within the reduced Hilbert space £d, 
where n+(x) = and n_(a;) = ±1 for each x. In the 
latter region {\ng\ ^ 1/4), on the other hand, the low- 
energy charge configuration should satisfy the condition 
ri+(a;) = for all x. Unlike the former case, the ground 
state of iJp is non-degenerate and forms a Mott insu- 
lator characterized by ni(x) — n2{x) = for all x. As 
Ej is turned on, the ground state of is mixed with 
the states with n^{x) — ±2. Accordingly, the relevant 
reduced Hilbert space is given by where n^{x) = 
and n^{x) — 0, ±2 for all x (see also Ref. Q). 

We first consider the region near the maximal- 
frustration lines, where we project the Hamiltonian 
Eq. (|l[) onto £d, and analyze the properties of the sys- 
tem near the maximal-frustration line (jrig — 1/2| <C 1/4), 
based on the resulting effective Hamiltonian. Given the 
projection operator P onto £d, the effective Hamiltonian 
up to the second order in Ej/Eq, 



H. 



QPM 



H, 



1 - P 



P, 



(5) 



can be obtained via the standard procedure . Imple- 
menting explicitly the projection procedure, we get the 
effective Hamiltonian describing a single spin-1/2 anti- 
ferromagnetic Heisenberg chain 

Heff=lJj2s%x)S^x+l) (6) 
- l^T. [S+{x)S-{x+l) + S-{x)S+{x+l)] , 

X 

where the exchange interaction and the uniaxial 
anisotropy factor are given by J = Ej/AEq and 7 = 
IQX^Ej /Ej, respectively. The pseudo-spin operators 
have been defined according to: 



S+{x) = Pe-'"^i(^^(l - P)e+*'^^(^'P 
S~{x) = Pe-"f''^''\l - P)e+*'^i(=")p. 

(Note the difference from the standard definition.) 



(7) 
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The effective Hamiltonian in Eq. (|^) includes contribu- 
tions from several complex processes, back in the charge 
picture: The first term in Eq. (||), which comes from the 
projection PH^P, simply describes the nearest-neighbor 
interaction of the charges in the form n-(x). On the 
other hand, the second term in Eq. (^, arising from the 
second order expansion PHj^^^HjP, describes the 
cotunneling process of two particles on different chains 
by way of a quantum mechanical virtual state with en- 
ergy of the order of Eq. Figure ^ shows schematically a 
particle at x on one chain and another particle at x+1 (or 
x—1) on the other chain hopping at the same time in the 
opposite direction. This cotunneling process plays the 
major role in the charge transport along the two chains, 
and drives the quantum phase transition in the system, 
as discussed below. 

The antiferromagnetic Heisenberg chain described by 
Eq. (^) has been extensively studied and is known, 
from the Bethe ansatz solution or the Sine-Gordon the- 
ory to exhibit a quantum phase transition at 7 = 1: 
For 7 > 1, it belongs to the universality class of the 
Ising chain in the renormalization group (RG) sense, 
and its ground state displays genuine long-range order 
in the staggered magnetization, i.e., {{—l)^S^{x)S^{0)) 
approaches a nonzero constant as a; — > 00. This long- 
range order in the staggered magnetization corresponds 
to the charge-density wave (CDW) in the charge picture 
of the original problem. For 7 < 1, on the other hand, 
the system described by Eq. is equivalent to the quan- 
tum XY chain, where the Mermin- Wagner theorem pro- 
hibits genuine long-range order. In this case, the system 
can be mapped to the repulsive Luttinger model, and 
the transverse component of the magnetization as well as 
the z-component of the staggered magnetization exhibits 



quasi- long-range order. Namely, both ((— l)^S'^(a;)S'^(0)) 
and (S'+(x)S'~(0)) decay algebraically with the distance 
X, and the system in the charge picture displays both 
the diagonal and the off-diagonal quasi-long-range or- 
der. This state may be regarded as the counterpart 
of the supersolid, possessing both the diagonal and off- 
diagonal (true) long-range order and proposed recently 
in 2D Josephson-junction arrays 

The properties of the repulsive Luttinger liquid phase, 
with both the diagonal and off-diagonal quasi-long-range 
order, has been discussed in Ref. [ pT| for a single 
Josephson-junction chain. It has been suggested that 



the system is extremely sensitive to impurities |15 and 
may make another insulator, different in nature from the 
CDW insulator. In the coupled chains, the repulsive Lut- 
tinger liquid phase has another remarkable feature of the 
current mirror. According to the basic transport mech- 
anism due to the cotunneling process of particles in the 
two chains, shown in Fig. ^ the current fed through one 
chain is accompanied by the secondary current in the 
other chain, with the same magnitude but in the oppo- 
site direction. Similar current mirror effects have also 
been pointed out for Ug = 0, where the mechanism is 
rather different and via the particle-hole pair transport 
(see below and Ref. ||^). 

We now turn to the region near the particle-hole sym- 
metry line. In the reduced Hilbert space £s with the 
condition n^(x) = satisfied, n_(x)/2 can be regarded 
as the number of particle-hole pairs located at x. The 
role of such particle-hole pairs can be analyzed by means 
of the imaginary-time path-integral representation of the 
partition function and its dual transformation Q. The 
Euclidean action, in the current-loop representation [Q, 
then reads 



S 



[nt{x,T) - Ug] \2CiCp^{x,x')\ [ni>{x\T) 



l.x.r 



(8) 



where the (imaginary) time has been rescaled in units 
of the inverse Josephson plasma frequency 07^7 ^ = 
h/^/4EjEj, the dimensionless coupling constant defined 
to be X = ^/EJ/16EI. Here {rn, Ji) may be viewed as 
the current in (1-1-1) -dimensions, satisfying the continu- 
ity equation 



Vrni{x, t) + VxJe{x, t) = 0. 



(9) 



With the capacitance matrix Eq. (^), it is convenient 
to decompose the action in Eq. into the sum S = 
S+ + S-: 



E, 



AKE 



-y 



n^lx, T) 



(10) 



where J±{x, t) = Ji {x, r) ± J2(x, t). 

The factor Eo / Ei in the component 5+ , which is enor- 
mous in the parameter regime of interest, again im- 
plies the condition nj^{x) — already mentioned. Fur- 
ther, the continuity equation in Eq. (^ requires J+{x) 
to be a constant on the average, which should obvi- 
ously be zero. Consequently, near the transition point, 
the system is effectively described by the action S- in 
Eq. (10), which is equivalent to the 2D XY model, and 
exhibits a BKT transition aX K = Kbkt ~ 2 /it [p|. 
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Here the transition, which is between the Mott insu- 
lating phase and the superconducting phase, is driven 
exclusively by the particle-hole pairs represented by the 
variable n_(x), whereas n+ and J+ merely renormalize 
the action 5_ and shift slightly the transition point [^. 
This shift of the transition point depends on the exter- 
nal charge Ug and may be estimated in the following 
way: The transition to the superconducting phase oc- 
curs when the Josephson-coupling energy Ej also over- 
comes the Coulomb blockade associated with a particle- 
hole pair. Since the Coulomb blockade increases with Ug, 
approximately given by SEorv^+iEi, the critical value of 
Ej is concluded to grow from the symmetry-line value 
IGK^j^rpE] as Ug is increased. It is also stressed that the 
BKT-type transition survives the gate voltage as long as 
the induced charge is sufRciently small {\ng \ <^ 1/4); this 
is in sharp contrast to the single-chain case, where break- 
ing the particle-hole symmetry by nonzero ng immedi- 
ately alters the universality class of the transition Q . 

In the parameter regimes other than those considered 
above, the behavior of the system may be inferred by 
the following argument: First, it is obvious that for 
Ej ^ Eq, the system should be a superconductor with 
each chain superconducting separately. Note that this 
superconducting phase, denoted by S, comes from the 
particle (Cooper pair) transport as usual, thus different 
in character from the superconducting phase in the re- 
gion Ej <^ Eq. In the latter, denoted by S', only the 
coupled chains as a whole is superconducting, with su- 
perconductivity arising from the particle-hole pair trans- 
port. Far away from both the particle- hole symmetry 
line and the maximal frustration lines, the Hamiltonian 
may be projected onto the subspace where ni{x) = 0, 1 
and n2{x) = 0,1 for all x, and single-particle processes 
dominate the transport in the system. 

The observations so far are summarized by the phase 
diagram displayed schematically in Fig. |[ The phase 
transitions of our main concern are represented by the 
thick solid lines, separating the CDW from the repul- 
sive Luttinger liquid (LL) and the Mott insulator (MI) 
from the superconductor (S'); the somewhat speculative 
boundaries discussed above are depicted by dashed lines. 
Here it is not clear within our approach whether the 
boundary between the repulsive Luttinger liquid region 
and the superconducting region in the phase diagram de- 
scribes a phase transition or merely a crossover. Further- 
more, even in the single-chain case, the properties of the 
repulsive Luttinger liquid phase is controversial, and the 
possibility of an intermediate normal phase has recently 
been raised as well [ p^ . 

Coupled chain systems can presumably be realized in 
experiment by current techniques, which have already 
made it possible to fabricate submicron metallic junction 
arrays with large inter-array capacitances [Q as well as 
large arrays of ultra-small Josephson junctions We 
also point out that quasiparticles have been safely dis- 



regarded in obtaining the equilibrium properties at zero 
temperature. 
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FIG. 1. Schematic diagram of the system. 
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FIG. 2. A typical second-order process via intermediate 
virtual states with energies of the order of Eq near the maxi- 
mal-frustration line. 
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FIG. 3. Schematic phase diagram of the coupled Joseph- 
son-junction chains. 
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